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Abstract 

We consider the large- U limit of the one-band Hubbard model at half-filling on a 
non-bipartite two-dimensional lattice. An external magnetic field can induce a three- 
spin chiral interaction at order 1/U 2 . We discuss situations in which, at low tem- 
peratures, the chiral term may have a larger effect than the Pauli coupling of electron 
spins to a magnetic field. We present a model which explicitly demonstrates this. The 
ground state is a singlet with a gap; hence the spin susceptibility is zero while the chiral 
susceptibility is finite and paramagnetic. 
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The ground state and low-energy excitations of Heisenberg antiferromagnets have 
been extensively studied in recent years. One reason for this interest is that the simplest 
model for strongly correlated electron systems, the one-band Hubbard model, reduces 
at half-filling to a Heisenberg antiferromagnet in the limit of large on-site Coulomb 
repulsion U. 1 Specifically, the Hubbard hamiltonian is 

H = - *y c t c 3 a + U n H t 1 ) 

17,(7 i 

where Uj = tj* is the hopping integral from site j to site i, a = ±1 or, 
equivalently, | and [ , and rii G = c\ G Ci a . At half-filling, J2i a Uicr 1S ec l ua l to 
the number of sites. In the large-?/ limit, the low-energy subspace of the total Hilbert 
space consists of states with exactly one electron at each site. To order l/U , this 
subspace is governed by the spin Hamiltonian 

Hsp = Jij ^2 (2) 

ij 

Here = (h/2) cJ CT t gg i c^i with the r being Pauli spin matrices (the indices 
a, a' are summed over). The couplings = 2 | tij | 2 jJJ are antiferromagnetic. 
At higher orders in l/U , the get renormalized and multi-spin interactions appear 
in H sp . These can be calculated as an expansion in tij/U by relating Eqs. (1) and 
(2) through an unitary transformation. 2 

In this paper, we examine what happens when this model is placed in a magnetic 
field . 3 The application of the field has two effects. Firstly, the Uj pick up a phase 
ex.p[(ie/ch) J? A-dr] where A is the vector potential. Hence the phase of the string 
Uj tjk • • • tki connecting three or more sites forming a closed curve is proportional 
to the flux enclosed by that curve. Secondly, there is the Pauli interaction given by 
v Si ■ ~B/h where v = — eh/ mc , with e (m) and c denoting the charge (mass) 
of the electron and the velocity of light respectively. (As explained below, the Pauli 
term takes the same form in the Hubbard and spin Hamiltonians). We are interested 
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in finding out whether the phases in tij induce any unusual terms in H sp , and what 
effects such terms may have. Since an external magnetic field breaks invariance under 
time-reversal T, we might expect H sp to reflect this. Namely, if B — > — B so that 
tij — >• t^* , there should be terms in H sp which reverse sign. It turns out that no 
T-violating terms are induced in H sp on bipartite lattices as one can show by using 
the particle-hole symmetry at half-filling. We transform Cj )(T — > a c\ _ a on the 
sites of one sublattice, and Ci, a — ► — a cj _ a on the other sublattice. This is a 
symmetry of (1) if tij — > t*j at the same time. Since Sj remains invariant under 
this transformation, H sp must be the same for B and — B . So the T— violation 
appears to lie entirely in the high energy subspace (states with one or more doubly 
occupied sites) for bipartite lattices. We also observe that on any lattice, particle-hole 
symmetry implies that H sp must remain invariant if tij — > — tji in (1). Hence if 
the t^ are all real, H sp cannot have odd powers of t^ . 

We must therefore consider non-bipartite lattices to obtain something interesting. 
The simplest example consists of three sites i,j and k forming a triangle. The 
perturbative expansion in tij/U is obtained by first writing the hopping term in (1) 
as the sum of three terms T , 7\ and T_i = , where T m increases the number 
of doubly occupied sites by m when it acts on a state. 2 The unitary operator relating 
(1) and (2) is then given by exp[ i K ] where K is a power-series in T m /U . (Note 
that the Pauli interaction commutes with all the T m and therefore with K). At half- 
filling, the low-energy subspace (states with no doubly occupied sites) is annihilated by 
both T and T_i since any hopping necessarily takes us to a state with one doubly 
occupied site. We then find that 

^^^SrB-i^^ + i T_i T T, (3) 

i 

When rewritten in the language of spin- 1/2 operators, the second term on the right 
hand side of (3) is the same as Eq. (2) while the third term is a three-spin chiral 
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interaction 3 of the form \x Si • Sj x Sk/h 3 where 

24 

\i = -^2 Im ( *v £ i fc £fci ) ( 4 ) 

This vanishes if the Uj are real. Let the magnetic flux enclosed by the triangle be 
$ . If we denote the magnitude | tjk tki I = ^ 3 for simplicity, then 

24£3 • f e$ ^ ^ 

One can estimate the relative magnitudes of this chiral term and the Pauli term 
for some 'typical' values of Uj, U and the area of the triangle A 4 For A = 2 A° 2 , 
the number eQ/ch is much smaller than 1 unless the field B reaches the fantastically 
large value of 10 4 Tesla. Hence we replace the sine in Eq. (5) by its argument, so 
that 

24t 3 eA n 

where 6 is the angle between B and the normal to the plane of the triangle. We then 
find that for t = 0.5 eV and U = 5 eV, the magnitude of the Pauli term is 
about forty times larger than the chiral term. This estimate follows from comparing 
the splitting in the ground state energy produced by the Pauli and chiral interactions 
for a triangle in which the three J^s are equal. Namely, we take the Hamiltonian on 
a triangle to be 

H S p — J [ Sj • Sj + Sj • Sfc + Sfc • ] 

(7) 

+ JL Si ■ Sj x S fc + \ ( Si + Sj + S fc ) • B 



If B = , the ground state of this Hamiltonian has a four-fold degeneracy with all 
states having total S = 1/2 . The magnetic field breaks this degeneracy completely 
with the Pauli and chiral terms contributing ±vB/2 and ±\/3/x/4 respectively. For 
the excited states with S = 3/2 , we observe that the chiral term has no effect. Thus 
the chiral interaction can only lower the energy of a state if it is no n- ferromagnetic. 



Although the Pauli term appears to be numerically much larger than the chiral 
term, one can think of two possible situations in which the chiral interaction dominates. 
The first example is one in which the ground state is a spin singlet and is chiral even in 
the absence of the magnetic field . We have in mind here the two-dimensional models 
discussed by Wen, Wilczek and Zee 5 where a spin Hamiltonian has two degenerate 
singlet ground states with opposite chiralities. We can say that each ground state has 
a non-zero chiral moment M c (defined later). Then an applied magnetic field picks 
out one of the two ground states due to an interaction of the form — M c B . The 
Pauli term S • B (where S = J2i ^ ) plays no role here because the ground states 
are singlets. However, these kinds of models often require a special choice of the two- 
spin couplings Jij as well as peculiar multi-spin interactions in order to produce the 
required ground state degeneracy. There are also papers which argue that frustrated 
antiferromagnets with only two-spin interactions can have chiral ground states. 5-7 
However, this has been questioned 8 and it seems to be quite difficult to have chiral 
ground states in the absence of an external magnetic field . 

The second example, which only has short-range two-spin interactions and does 
not require a fine-tuning of the couplings, is one in which the ground state is a singlet, 
unique and non-chiral in the absence of the magnetic field . Further, there is a gap A to 
states with total spin greater than zero. Then the ground state continues to be a singlet 
in the presence of a field if | vB | <C A. But it may develop a chiral moment M c to 
first-order in B and one can define a chiral susceptibility Xc = ( dM c /dB )b = o ■ So 
one may have a finite Xc even though the spin susceptibility Xs = 9 { S ) / dB = . 
We now present a two-dimensional model which explicitly demonstrates all this. Inci- 
dentally, it is the only two-dimensional spin model that we are aware of in which the 
ground state and the low-lying excitations can be found exactly (for the Hamiltonian 
Hq given below). 

Our model, shown in Fig. 1, consists of chains of rhombuses which are coupled to 
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each other in the form of a brick lattice. Each rhombus is formed out of two triangles 
with a common base. The number of rhombuses is N/3 if the number of sites is 
N . Starting from a Hubbard model with only nearest-neighbor hoppings, the spin 
Hamiltonian in a magnetic field is given up to order l/U 2 by 

C 

H sp — Ho + \i —3 
where H = J v + \ J2 S * " B (8) 

ij ' ■ i 

and C = S 3 • S 4 x ( S 2 - S 5 ) 

a 

The chiral term C is a sum over rhombuses labelled by the index a , with each rhombus 
contributing the sum of two three-spin terms as indicated in Eq. (8). (See Fig. 1 for 
the site labels 1 to 6 in and around a typical rhombus). Note that the Pauli term 
has been included in an 'unperturbed' Hamiltonian H , while the chiral term will be 
considered perturbatively in the following. If we define parity to be the transformation 
which exchanges the top and bottom sites of the vertical bonds inside all the rhombuses 
simultaneously (namely, 3 <-> 4 in the figure), then C is odd under parity while H 
is even under parity. 

In H sp , the couplings Jij on the vertical bonds inside the rhombuses, the slanted 
bonds, and the vertical bonds joining the chains are denoted by J\ , J2 and J3 
respectively as shown in Fig. 1. Let us assume that J\ > 2J2 and that | vB \ is much 
less than both ( J\ — 2J 2 )h and J3 Ti . Then one can prove that the ground state 
of Hq is a singlet, unique and has a gap Ao to all excitations. Let us first introduce 
the notation Oij and 1^ for the singlet and triplet states respectively formed from 
the spins at sites i and j. Note that Oij = l/\/2 ( | i | j [) — \ i [ it}) is 
antisymmetric under an exchange of i and j, while the three states collectively denoted 
by lij are all symmetric. Then the ground state of Hq is the state ^0 given 
by the product of singlets • • • O12 ® O34 ® O56 • •• following the labels in Fig. 
1. Namely, each of the vertical bonds form a singlet. The ground state energy is 
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E = - ( 7V/8 ) ( 2 Jx + J 3 ) . 

To prove that ipo is the ground state and that there is a gap A , let us 
write Hq = Hi + H 2 where Hi is the same as H except that the cou- 
plings J\ are replaced by J\ — 2J 2 and the couplings J 2 are replaced by zero. 
Thus Hi is a sum of disconnected two-spin Hamiltonians involving only the verti- 
cal bonds, while H 2 is a sum of disconnected four-spin Hamiltonians of the form 
( J 2 /2) [ ( S 2 + S 3 + S 4 ) 2 + ( S 3 + S 4 + S 5 ) 2 ] for each rhombus. It is then easy 
to find the complete spectra for both Hi and H 2 . For Hi , ipo is the unique ground 
state and there is a finite gap Ai = min ( Ji — 2 J 2 — | vB | , J3 — | vB | ) to the 
space of states orthogonal to ipo • For H 2 , ipo is a ground state but there is no gap 
to its orthogonal subspace. It then follows that ipo is the ground state of Hi + H 2 
and there is a gap A > Ai to all other states. 9 We also note that the total spin 
on any bond of type Ji (e.g. (S3 + S4 ) 2 ) commutes with Ho , so that there are 
N/3 operators which can be diagonalized along with Hq . This important property 
proves to be very useful. For instance, it implies that a low-lying excitation can only 
have a finite number p of Ji bonds forming triplets. (Such a state is separated from 
the ground state by a gap A > p ( Ji — 2J 2 — \ vB | ) by a similar argument 
involving Hi and H 2 ). Further, such an excitation can only differ from the ground 
state in a local neighborhood of those p triplet bonds or due to some isolated J3 
bonds forming triplets instead of singlets. Hence, all low-lying excitations are localized 
and dispersionless (with energy independent of the momentum). 

One can check that -00 is non-chiral, namely, ( ipo \ C \ ipo ) = . The simplest 
way to show this is to write C = C a , where C a is the sum of the two chiral 

terms in a rhombus a. Consider the rhombus made out of sites 2, 3, 4 and 5 in Fig. 
1. Both ipo and C a for that rhombus are odd under the parity transformation 
3^4. Hence C a \ ipo ) is even under parity, i.e., it contains I34 rather than 
O34 . Hence ( -00 I C a \ ipo ) must be zero. We will now assume that the term 
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[iC in (8) only changes the ground state perturbatively because we expect the gap to 
survive for a finite range of ji around ji = . We can therefore use second-order 
perturbation theory to compute the ground state energy E (fi) to order /x 2 . Since fi 
is proportional to B , this will give us the chiral moment M c = — dE Q ((i)/dB and 
the chiral susceptibility \c = dM c /dB . 
The second-order expression is 

E ( P ) - E a = ^ £ £ (C ° )o »_ (C % U (9) 

n^O a,p 

where (C a ) mn = { tjj m \ C a \ ip n ) and m, n label the eigenstates of H . Now, we 
know that the state Cp \ tpo ) has the Ji bond in rhombus (3 forming a spin triplet 
while the J\ bonds in all other rhombuses are singlets. This implies that all the terms 
in (9) with a ^ (3 must vanish. Next let us consider a particular rhombus labelled 
P and the six sites in and around that rhombus as labelled in Fig. 1. Then (Cp ) n o 
can only be non-zero for a finite number of states ip n since the singlet subspace of 
those six spins is five-dimensional. Also, (Cp ) n o can be non-zero only if the vertical 
bond in the rhombus (3 is a triplet in the state tjj n . An explicit calculation shows 
that (Op ) n o is actually non-zero for only two states for which E n — E are given 

by 

E± - E = Ji - J 2 + 5 J 3 ± ^ J 2 2 + ^ (10) 

respectively. ( E± and E are independent of the magnetic field B since all the states 
being considered are spin singlets.) We eventually find that 

T-i / \ T-i N 2 Ji - J 2 + 2J 3 

EM - E ° = -J" (E + -E )(E_-E ) M 
This is of order t 4 /U 3 since the couplings Jij ~ t 2 /U . This and Eq. (6) yield a 
paramagnetic susceptibility 

N ( 24t 3 eA \ 2 Ji - J 2 + 2J 3 

Xc = 4 { JP Hi COS9 J (E + - E )(E_ - Eo) (12) 
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to order t 4 /U 3 . 

We emphasize that (11) is not the complete expression for the ground state energy 
to order 1/U 3 , since we have ignored terms of order 1/U 3 in deriving H sp in (8). 
These terms do contribute to the energy in first-order perturbation theory. However, 
one can show that these terms are not chiral because they are even under parity. They 
are independent of \i and hence do not contribute to the chiral quantities M c and Xc 
to order t A /U 3 . 

At finite temperatures, this model no longer has Xs = . But if the temperature 
is small compared to the gap, then Xc will continue to be much larger than Xs ■ 

Our model is somewhat peculiar in that the two-spin correlation in the ground 
state of H is exactly zero beyond a short distance. However, this property is unlikely 
to survive once we take into account the chiral terms and higher order terms in 1/U 
which couple spins on no n- neighboring sites. We then expect the two-spin correlation 
to go to zero exponentially at large separations because of the gap above the ground 
state . The ground state is therefore a spin-liquid which is dominated by short-range 
valence bonds. 

More realistic models which do not have a gap to spin excitations will generally 
have both Xc and Xs non-zero even at zero temperature. For instance, one can 
consider a Hubbard model on a triangular lattice, or on a square lattice with both 
nearest-neighbor and next-nearest-neighbor hoppings. Whether Xc will be comparable 
to or much smaller than Xs will then depend on the properties of low-energy excitations 
in the absence of the magnetic field . For instance, if there are singlet chiral states 
lying very close to a non-chiral ground state, then one would expect Xc to be large. 
An important (and perhaps experimentally observable) difference between the two 
susceptibilites is that Xc depends on the orientation of the magnetic field with respect 
to the plane containing the sites of the spins. 

To conclude, we have seen that a spin system which arises from an underlying 
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Hubbard model can develop chiral interactions when placed in a magnetic field. Al- 
though these interactions are small, they may lead to an interesting low-temperature 
phase resembling a chiral spin-liquid. 
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Figure Caption 

1. The model showing six sites labelled 1 to 6 in and around a typical rhombus. The 
three different antiferromagnetic couplings J\ , J2 and J3 are also indicated. 
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